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ABSTRACT 

Recent determinations of the upper mass limit to the local initial mass function (IMF) 
claim a value of my = 50 ±10 M , based upon direct comparisons of the observed oxygen 
and iron abundances in metal-poor stars with the predicted stellar yields from Type II 
supernovae (SNe). An unappreciated uncertainty in these analyses is the input physics 
intrinsic to each SNe grid, and its effect upon stellar nucleosynthesis. We demonstrate 
how such uncertainties, coupled with the uncertain metal-poor halo star normalization, 
while allowing us to set a lower bound to my of ~ 40 M , nullifies any attempt at 
constraining the upper bound. 

Subject headings: nuclear reactions, nucleosynthesis, abundances — stars: luminosity, 
mass function — stars: supernovae — galaxy: evolution 



1. Introduction 

The initial mass function (IMF), a measure of the 
distribution of masses, at formation, of a given stellar 
generation, is one of the key components of galaxy 
evolution modeling. Despite its importance, accurate 
determination of the IMF remains one of the most 
elusive problems in modern astronomy. This elusive- 
ness is manifested both in ongoing attempts to under- 
stand its physical underpinings (e.g. Padoan, Nord- 
lund & Jones 1997), as well as simply characterizing 
its shape and mass limits from an observational tack 
(e.g. Kroupa, Tout & Gilmore 1993; Scalo 1986). 

In its simplest form, the IMF can be considered 
a power law of the form rruj)(m)dm oc m _ ' 1+I 'dm, 
where m(j)(m)dm is the number of stars born in the 
mass interval m — ► m + dm. The goal for theorists 
and observers, alike, is then the determination of the 
slope of this function well as its upper and lower 
limits (mu and m e , respectively). 

While the slope of the IMF, at least in the solar 
neighborhood (and for masses greater than a few solar 
masses), would appear to lie somewhere in the range 
x w 1.3 (Salpeter 1955) to x ~ 1.7 (Kroupa, Tout 
& Gilmore 1993), and the lower mass limit is close 
to m e w 0.2 M (Bahcall et al. 1994), the upper 
mass limit totj still remains highly uncertain. Even 
a cursory examination of the literature corroborates 
this point, with values in the range totj w 20 — > 200 
M suggested by a variety of direct and indirect 
techniques (e.g. Maeder & Meynet 1989; Klapp & 
Corona-Galinda 1990; Pagelet al. 1992; Maeder 1992; 
Massey et al. 1995; Kudritzki 1997). 

A hybrid approach to determining totj, combining 
predictions of the theoretical yields from Type II su- 
pernovae (SNe), with the observed abundances in the 
old metal-poor stars (i.e. those which bear the clear 
imprint of yield "pollution" from these same SNe, 
with no "dilution" from Type la SNe, whose progeni- 
tor lifetimes are considerably longer than the Type II 
timescales), has been the subject a recent series of pa- 
pers (Tsujimoto et al. 1995,1997; Yoshii, Tsujimoto 
& Nomoto 1996). The premise here is that because 
Type II SNe a-element (e.g. O, Mg, Ne) yields are 
a strong function of progenitor mass, whereas prod- 
ucts of explosive burning (e.g. Fe, Si, Ca) are less so, 
the IMF- weighted average of their ratios must neces- 
sarily also depend strongly upon x and m\j. Tying 
these yield "averages" to the halo abundances then, 
in principle, provides a unique indirect probe on the 



upper mass limit to the IMF.Q 

Following this technique, Tsujimoto et al. (1997) 
recently concluded that the upper mass limit to the 
IMF in the solar neighborhood is mu = 50 ± 10 Mq . 
To do so, they made explicit use of the Tsujimoto 
et al. (1995) compilation of Type II SNe yields. What 
was not fully appreciated in their study though was 
just how dependent their result was to this particular 
yield compilation and the adopted halo abundance 
normalization. It is to this lack of appreciation that 
our current study is addressed. 

After providing a minimal introduction to the 



model ingredients in Section 2.1, we demonstrate in 



Sections 2.2 and 2.3 that this technique results in 
totj = 50 ± 10 Mq only for the Tsujimoto et al. 
(1995) yields, combined with a halo normalization 
of [0/Fe]h = +0.41. Duplicating the analysis with 
"competing" yield compilations which sample a wide 
variety of convection and mass-loss treatments (e.g. 
Woosley & Weaver 1995; Langer & Henkel 1995; Ar- 
nett 1996), clearly demonstrates that Tsujimoto et al. 
(1997) have significantly underestimated the uncer- 
tainty associated with thir determination of mj . Our 
results are summarized in Section |^. 

2. Analysis 

2.1. The Basic Formalism 

The observed oxygen-to-iron abundance ratio in 
halo dwarfs and giants, [0/Fe]h, can be linked (see 
Tsujimoto et al. 1997 - hereafter, T97 - for details) 
to the theoretical Type II SNe yields m c A and m^ e via 



xl0 [o/Fc ]h = /io" m o m 

/™ u m^m-(i+-)dm' 
(1) 

where the coefficient (0/Fe) Q =7.55 is the solar (me- 
teoritic) mass fraction ratio from Anders & Grevesse 
(1989), and x is the slope of the IMF (recall Section 

0)- 

By specifying the low metallicity "plateau" value 
for [0/Fe]h (e.g. Figure 11 of Timmes, Woosley & 
Weaver 1995; Figure 2 of Bessell, Sutherland & Ruan 



1 Tsujimoto et al. (1997) extend this "hybrid" approach to si- 
multaneously constrain the IMF mass limits my and mp, as 
well as the slope x. We shall only be concerned with the m\j 
determination in what follows, primarily for brevity, but also 
because the lower mass constraint rest squarely upon uncertain 
photometric calibrations. 



2 



1991), and the IMF slope x, equation [l] provides a 
unique upper mass limit my, for a given Type II SNe 
yield compilation. 

2.1.1. Ingredients 

Essentially, there are only two input ingredients 
which interest us, in what follows. First, the source 
of Type II SNe yields (right-hand side of equation [l]) , 
and second, the halo's "plateau" [0/Fe]h (left-hand 
side of equation [l]) . Let us briefly comment on the 
latter ingredient first. 

There is no debate as to the reality of the a- 
element overabundance seen in halo dwarfs and giants 
(i.e. [0/Fe]h <; +0.2), but there does remain a fac- 
tor of ~ 2 uncertainty as to its asymptotic "plateau" 
value for [Fe/H]< -3. T97 adopt [0/Fe] h = +0.410 
but it is apparent that this depends somewhat upon 
the sample selection. Visual inspection of Figure 
11 of Timmes et al. (1995), Figure 2 of Bessell 
et al. (1991), and Figure 1 of T97, each drawn from 
slightly different sources, shows that the plateau lies 
roughly between [0/Fe]« +0.4 -> +0.6; i.e. T97's 
value of +0.41 appears to lie on the lower end of 
the plateau "distribution" . T97 did not demonstrate 
how the predicted mjj changes as a function of this 
halo "normalization", a point to which we return in 



Section 2.3, where we duplicate their analysis using 
[0/Fe]h = +0.60 (i.e. that favored by Bessell et al. 
1991). 

The Type II SNe yields are the primary input into 
equation |l|. T97 adopt their earlier 1995 yield compi- 
lation (T95, hereafter), which itself is an offshoot of 
the Thielemann, Nomoto & Hasimoto (1996) models. 
The solid curves in Figure ^ represent the T95 oxygen 
and iron yields adopted in the T97 analysis. 

It is readily apparent though, as even a cursory 



2 T97's claim to adopt [O/Fejh = +0.41 is not entirely correct, 
as their code does not follow O and Fe, per se, but only the 
isotopes 16 O and 56 Fe, thereby underestimating the true iron 
yield by ~ 10% (and oxygen by a smaller amount). This has the 
effect that instead of using (O/Fe)0 = 7.55 in equation hi, T97 
use (O/Fe) = 8.20 (i.e. ( 16 O/ 56 Fe) = 8.20). By itself, this 
is not a problem; the problem arises in that T97 then normalize 
equation FJ to the observed halo ratio [0/Fe] h = +0.41 - this 
observed ratio, though, reflects the total O and Fe abundances, 
and not just the isotopes ie O and 56 Fe. In practice, what 
this means is that T97 use (O/Fe)^ = 21.1 (recall equation 
E), which corresponds to [0/Fe] h = +0.45, and not +0.41. To 
truly normalize their models to [0/Fe]h = +0.41, T97 should 
have used (O/Fe)^ = 19.4 and included the total O and Fe for 
the yields in the right-hand side of equation hi 
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Fig. 1. — Oxygen (upper four curves) and iron (lower 
three curves) yields as a function of progenitor mass, 
for the primary yield compilations considered here 
- T95=Tsujimoto et al. (1995); LH95=Langer & 
Henkel (1995); A96=Arnett (1996); WW95=Woosley 
& Weaver (1995); M92=Maeder (1992). The models 
of M92, LH95 and A96 are only evolved up to the 
completion of oxygen burning, hence an iron yield is 
not directly associated with either grid. We indirectly 
associate Arnett's (1991) iron yields with his newer 
1996 oxygen yields. See text for details. 



glance at Figure 1 should show, that there still remain 
substantial uncertainties in the yield computations for 
Type II SNe. These differences have already been de- 
scribed in detail by Arnett (1995), Langer (1997), and 
Gibson, Loewenstein & Mushotzky (1997), to name 
just a few. 

T97 assume that these differences can be rec- 
onciled by considering what they (mistakenly) as- 
sume is the maximal deviation from their canonical 
(T95) yields. This deviation is taken to be a mass- 
independent 30% excess added to their oxygen yields; 
iron, though, is assumed to be invariant. T97 jus- 
tify this 30% excess as the "extreme" , since Woosley 
& Weaver's (1995) (WW95, hereafter) [Fe/H]=+0.0 
oxygen yields are roughly 30% greater than those of 
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T95 (also [Fe/H]=+O.O).0 Herein lies several (inter- 
twined) problems. 

First, solar metallicity yields are inappropriate for 
the discussion at hand; T97 attempt to constrain 
mjj by comparing the abundance ratios in metal-poor 
stars (i.e. [Fe/H]< —2) with the yields predicted from 
[Fe/H]=+0.0 Type II SNe models. A more appro- 
priate course of action is to consider a much lower 
metallicity grid of models. Because T95 do not com- 
pute [Fe/H]<+0.0 models, T97 could not explicitly 
consider metallicity effects. 

If we compare (properly) T95's oxygen yields with 
WW95's sub-solar metallicity grid, the 30% "excess" 
adopted by T97 is no longer relevant. Figure [l] shows 
that WW95's [Fe/H]=-4.0 models have oxygen yields 
virtually indistinguishable to the [Fe/H]=+0.0 mod- 
els of T95, except for m > 35 M , where WW95 now 
lies systematically below T95. 

These subtle differences in the WW95 and T95 
oxygen yields are interesting in their own right (e.g. 
Langer 1997), and are certainly encapsulated by the 
30% error budget adopted by T97. Unfortunately, a 
more fundamental flaw made by T97 was in assum- 
ing that all of the uncertainty in the oxygen yields 
could be accounted for by this 30% factor. The prob- 
lem with assuming that the T95 oxygen yields form a 
lower envelope, with WW95 forming the upper, can 
best be appreciated by referring to the two other pri- 
mary sources of SNe models, besides those of T95 and 
WW95 - i.e. Langer & Henkel (1995, hereafter LH95) 
and Arnett (1996, hereafter A96). 

In Figure [ll we can see that both LH95 and A96 
agree (roughly) with one another, ^ but lie a factor 
of - 3 below the T95 and WW95 predictions. This 
factor of three uncertainty in the oxygen yields is a 
far cry from the 30% uncertainty advocated by T97, 
but is entirely in keeping with that found by Langer 
(1997). This latter reference is an invaluable resource 
for those interested in properly assessing the origin 
of the large uncertainties in different modelers' oxy- 
gen yields. In the same vein, a useful comparison 
of the primary differences in the model input physics 

3 Actually, WW95's [Fe/H]=+0.0 oxygen yields (their favored 
Case B energetics) are greater than those of T95 only for m < 
30 Mq . 

4 The treatment of convection, the adopted 12 C(a,7) 16 re- 
action rate, and the initial evolutionary state of the models, 
were all significantly different in T95, than in either LH95 or 
A96. LH95 also considered sub-solar metallicity models and 
self-consistent mass-loss, both of which the others did not. 



can be found in Table 2 of Gibson, Loewenstein & 
Mushotzky (1997). At this point, we are not advo- 
cating one grid of oxygen yields over another, simply 
drawing attention to the fact that T97 have underes- 
timated its uncertainty. 

Figure [l] reflects the present-day state-of-the-art, 
as far as oxygen yields goes, and can be consid- 
ered to be the modern equivalent of Wang & Silk's 
(1993) Figure 1. The similarity of the Arnett (1978), 
Woosley & Weaver (1986), and Thielemann, Nomoto 
& Hashimoto (1994) oxygen yields, reflected in the 
latter figure, is not in dispute, but it should also be 
readily apparent (from our Figure [l]) that supplement- 
ing the current analogs of these older models (i.e. 
A96, WW95, and T95, respectively) with the new 
models of LH95 demonstrates that the agreement is 
no better than a factor of ~ 2 — > 3 at a given initial 
mass, in agreement with that found by Langer (1997). 

Figure || parallels that of Figure [j} but restricts it- 
self to the oxygen yields from the published compila- 
tions of Arnett (1978,1991,1996), hereafter A78, A91, 
and A96, respectively. The subscript 'He' or 'ZAMS' 
denotes the initial evolutionary state of the stellar 
models - pure Helium core or zero age main sequence. 
The former implies that a ZAMS mass-He core mass 
relation has been applied a posteriori, similar to the 
models of T95, whereas the latter avoid such com- 
plications by modeling the evolution self-consistently. 
We include a parallel analysis using each of the Ar- 
nett generations as (i) A78 has already played a large 
part in the analysis of Wang & Silk (1991), and for 
completeness sake should therefore be included here; 
(ii) A78 and A91 are identical in all respects, save the 
adopted He core mass-ZAMS mass relation, a useful 
comparison in and of itself; (iii) A96h c and A96zams 
are similar in all respects, save initial evolutionary 
state, thereby allowing us to better appreciate this 
crucial difference in the modeling. Langer (1997) has 
previously alluded to the danger of comparing yields 
from models evolved from the ZAMS versus those 
evolved from He cores; we are now in a position to 
quantify this "danger" , as applied to the determina- 
tion of the IMF upper mass limit. 

A further assumption inherent to the T97 analysis 
is that the iron yields are invariant and best repre- 
sented by those of T95. Iron predictions are particu- 
larly problematic as the exact location of the mass-cut 
leads to enormous uncertainties (Thielemann et al. 
1996). For the sake of self-consistency in the input 
physics though, we have not made the same a priori 
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Fig. 2. — Oxygen yields as a function of progeni- 
tor mass, from the four compilations of Arnett (1978, 
1991, 1996). Only the A96zams grid was evolved 
self-consistently from the zero age main sequence 
(ZAMS), while the remainder were evolved helium 
cores with an assumed a posteriori ZAMS mass-He 
core mass relation. 

assumption regarding iron's invariance, and have at- 
tempted to match the oxygen predictions with iron 
from the same grid. This is easy in the case of T95 
and WW95, as they follow the evolution right through 
core-collapse and explosive nucleosynthesis. Of the 
different Arnett grids described herein, only A91 in- 
cluded a specific entry for iron, so that is used for 
each of A78, A91, and A96, hereafter. Proper self- 
consistency with Arnett 's yields will be restored after 
he applies core collapse and explosive nucleosynthesis 
to the new A96 model grid. 

We are now in a pos itio n to anticipate some of the 
conclusions of Section |2.2| . To do so, we weight the 
mass-dependent oxygen and iron yields of Figures [j] 
and H by some canonical IMF, which we will take 
to be of the slope favored by Salpeter (1955) - i.e. 
x = 1.35 - with Type II SNe progenitors assumed 
to span the mass range 10 — > 50 M . The resulting 
mean IMF-weighted masses of oxygen and iron (and 
their logarithmic ratios relative to the solar ratio) - 
i.e. < O >, < Fe >, and < [O/Fe] > - for each of the 



yield "pairs" discussed thus far, are listed in Table 
[I]. The subscripts ZAMS and He refer to the initial 
evolutionary state of the models - i.e. zero age main 
sequence or simple helium cores. 

The factor of ~ 3 uncertainty due to oxygen (and 
to a lesser extent, that due to iron) is immediately ap- 
parent in columns 3 and 4, their combination leading 
to a factor of ~ 2 uncertainty in the IMF-weighted 
[O/Fe]. 

Recalling that T97, using the T95 yields, assumed 
a metal-poor normalization for equation |l] of [O /Fe] h « 
+0.4, and subsequently found m\j ~ 50 M , it should 
not be surprising to see that our entry in Table |l| for 
T95, generated using mu — 50 M Q , has an IMF- 
weighted [O/Fe] similar to T97's adopted halo value. 
Conversely, the A96 entries lie substantially below 
this halo value. Because [O/Fe] increases with in- 
creasing progenitor mass (see Figure |l|), it should be 
apparent that to recover an [0/Fe]« +0.4, with the 
same IMF slope, will require a greater m\j. 

2.2. Quantifying m^'s Dependence Upon Yield 
Source 

For a given halo "plateau" normalization [0/Fe]h, 
the vari ous yi eld compilations of Figures and I (and 
Section 2.1.1 ) can now be used in conjunction with 
equation |l| to determine how the upper mass limit 
my varies as a function of IMF slope x. The results, 
for a variety of models, are shown in Figure |^. Let us 
restrict ourselves, for the time being, to the six curves 
which were generated assuming the default halo nor- 
malization of [O/Fe] =+0.41. 

The first model to note from Figure || is that given 
by the solid curve, which effectively recovers the re- 
sults of T970 - i.e. mu w 48 -> 62 M , for reason- 
able IMF power-law slopes spanning x ~ 1.2> (Salpeter 
1955) to x ss 1.7 (Kroupa et al. 1993). Ag ain, r ecall- 
ing our comments at the end of Section 2.1.l| , this 
low predicted value for m\j was already anticipated, 
based upon the T95 yields possessing the greatest 
IMF- weighted <[0/Fe]>, for the "standard" IMF em- 
ployed in constructing Table []]. 

The most important result to take from Figure ^, 
though, is not so much the behavior of the T95 curve 
alone, but the parallel behavior for the models em- 
ploying the alternate Type II yields - i.e. those de- 



5 The footnote to Section [z.l.l] clarifies why there 
residual ~ 10% offset between the two. 



remains a 
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Table 1 
IMF-Weighted Yields 3 - 



Source* 3 


[Fe/H] c 


< o > d 


< Fo > cl 


< [O/Fe] > 


T95 Hc 


+0.0 


1.810 


0.091 


+0.42 


M92 ZA MS 


-1.3 


1.508 


n/a 


n/a 


WW95zams 


-4.0 


1.445 


0.087 


+0.34 


LH95zams 


-1.0 


0.841 


n/a 


n/a 


A78° +A9lj* 
A9l8+A9l| 
A96» c +A9l| 
A 96g AMS +A9l| c c 


+0.0 


1.711 


0.068 c 


+0.52 


+0.0 


1.189 


0.068 6 


+0.37 


+0.0 


0.983 


0.068 6 


+0.28 


+0.0 


0.678 


0.068 6 


+0.12 



a IMF slope x = 1.35, over the range 10 -> 50 M Q . 
b Source of oxygen yields. 

c Metallicity of models from which oxygen yields were derived. 

d IMF-weighted yield mass, in Mq. 

e Arnett (1991) [Fe/H]=+0.0 iron yields. 



noted WW95 and A960 in the legend to the figure. 
For a Salpeter (1955) IMF slope (i.e. x = 1.35), the 
predicted upper mass limit for the T95 yields was 
totj « 49 Mq . For comparison, the corresponding totj 
for the WW95, A78 and A96 Z ams yields is 61, 35, and 
131 M , respectively, with the A91 and A96h c values 
lying between A78 and A96zams- Because the depen- 
dence of totj upon IMF slope x is stronger in the non- 
T95 yields (especially for A96zams), for the Kroupa 
et al. (1993) slope (i.e. x = 1.7) mu is more strongly 
affected. The value favored by the T95 compilation 
(i.e. 62 M©) is replaced, in this case, by values rang- 
ing from - 95 M , for WW95,[] to > 200 M , for 
A96zams- 



6 We have not shown the results for the LH95 oxygen yields in 
Figure |i| as they parallel closely those of A96zams ■ 

7 Because the maximum mass considered by WW95 was 40 Mq, 
we make the (perhaps mistaken) assumption that linear extrap- 
olation beyond the grid extreme holds. For example, such an 
assumption leads to a predicted oxygen yield of 14 Mq for an 
m = 100 Mq star. If we were to arbitrarily double this yield 
to ~30 Mq, and set the iron yield to a constant 0.075 Mq for 
m > 35 Mq, we would find an approximate 5%/15% reduction 
in the predicted mrj, for IMF slopes x = 1.3/1.6. Even more 
extreme, if we increase the same oxygen yield to 50 Mq, and 
set the iron yield to zero for m > 35 Mq, we would find an ap- 
proximate 15%/30% reduction in the predicted mu, for IMF 
slopes x = 1.3/1.6. Obviously, self-consistent yield predictions 
for the WW95 grid, for masses m > 40 Mq, are required, al- 
though the previous examples do provide a feel for maximum 
effect involved. 



To summarize, only if one assumes that the T95 
yields are the definitive representation of low-metallicity 
massive star nucleo synthetic yields can one conclude 
that the upper mass limit to the solar neighborhood 
IMF is totj = 40 — > 60 Mq. By adopting, in turn, 
the three other primary yield sources (i.e. WW95, 
LH95, and A96), we sample a much fairer represen- 
tative input physics "parameter space" , leading to 
the more conservative result that we are able to con- 
strain the lower bound to the upper mass limit to be 
~ 40 M Q , but that the upper bound remains effec- 
tively unconstrained. If we force the Salpeter (1955) 
IMF slope (i.e. x = 1.35) to hold, the uncertainty 
in the yields still only allows us to constrain totj to 
the range ~ 40 — > 140 M . This yield- dependence of 
totj 's determination was not appreciated in the origi- 
nal analysis by T97. 

It is interesting to question how one might recover 
a result of totj k 40 — » 60 M , when adopting the 
A96zams yields. Assuming the oxygen as inviolate, 
it is apparent from inspection of Table |l| that we 
are required to reduce the IMF-weighted iron yield 
by a factor of ~ 2. This cannot be of the form of 
a blanket, mass-independent reduction, since we are 
constrained by the observations of SNe 1987a and 
1993j (Thielemann et al. 1996). Only by setting 
the iron yield to zero for to < 12 M and to > 30 
M (i.e. retaining the 15, 20, and 25 M iron pre- 
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Fig. 3. — Upper mass limit m\j, as a function of 
IMF slope x, for the Type II SNe yield compilations 
of Figures 1 and 2. 

dictions of A91, but ignoring the contribution from 
other masses), can we reduce the A96zams predic- 
tion for mu from ~ 130 — * 200 Mq to that found 
using T95 (i.e. m v « 40 -► 60 M Q ). 

In passing, we note that linear interpolation for 
both oxygen and iron were assumed in equation [l]. 
If one adopts a logarithmic interpolation, as Yoshii 
et al. (1996) did, the predicted values of my would 
be reduced by ~ 10 — ► 25%, from those shown in 
Figure |. 

2.3. Quantifying rou's Dependence Upon the 
Halo Normalization 

T97 did not consider how their conclusions might 
be influenced by the adopted mean "halo" [0/Fe]h 
(equation [l]) . Recall that the solid curve in Figure || 
supported my ~ 48 — > 62 M Q , for IMF slopes x « 
1.3 — » 1.7, for the T95 yields and a normalization 
[0/Fc] h = +0.41. 

If instead of adopting T97's favored halo normal- 
ization, we took that from Bessell et al. (1991) - i.e. 
[0/Fe] h = +0.60 - we would recover the "T95;+0.60" 
curve of Figure 0. In other words, this ~ 0.2 dex (i.e. 
~ 60%) higher normalization, for the T95 yields, in- 
creases the predicted mu by ~ 60% (from 48 M to 



75 M , for x = 1.3) to ~ 140% (from 62 M to 148 
Mq, for x = 1.7). The same conclusion is reached 
when we adopt the WW95 yields, as shown by the 
"WW95;+0.60" curve of Figure || A general rule-of- 
thumb for Salpeter (1955) IMF slopes (i.e. x = 1.35) 
is that a given percentage increase in the halo nor- 
malization, is accompanied by the same percentage in- 
crease in the predicted mu, regardless of yield source. 
For slopes x ~ 1.7, the increase in mu with increasing 
normalization is generally ~ 2 — > 3x greater. 

If we simply take the halo normalization to be 
[0/Fe] h = +0.5 ± 0.1, and accept T97's allowable 
range of IMF slopes (i.e. x w 1.3 — > 1.6), Figure 
^ would lead us to conclude that mu ~ 40 — ► 140 
Mq better represents the valid range of mu, for the 
T95 yields- we feel that this would have been a more 
realistic range than the mu ~ 40 — > 60 Mq claimed 
by T97. 

3. Summary 

T97 have recently revitalized interest in using IMF- 
weighted Type II SNe yields as a direct probe of 
said IMF's upper mass limit mu, by comparison with 
the observed abundance ratios in metal-poor Galactic 
stars. The beauty of this technique lies partially in its 
simplicity - for a given IMF slope, there is effectively 
only one free parameter - the yield source. Adopting 
the T95 yields, T97 found m v = 40 -> 60 Mq, for 
reasonable IMF slopes (i.e. x = 1.3 — > 1.6). 

The primary concern we have regarding T97's anal- 
ysis lies in their underlying assumption that all of the 
uncertainty in the stellar model physics could be en- 
capsulated in a 30% (0%) error budget for oxygen 
(iron). It should be obvious from Figures |l| and |[ 
and Langer (1997), for example, that this assumption 
is incorrect, and that a more realistic error budget 
would allow for up to an order of magnitude greater 
leeway. Stellar models are simply not developed to 
the level that is inherently assumed by T97 - convec- 
tion, overshooting, mass-loss, reaction rates, metal- 
licity, C/O-core masses, fallback onto the remnant - 
each conspire to increase the uncertainties to the de- 
gree reflected by the yields shown in Figures | and 

I 

A secondary concern is T97's inherent assumption 
that the halo normalization [0/Fe]h = +0.41 has 
no associated uncertainty. Since values as high as 
[O/Fe]=+0.6 are still favored by some, this ~ 60% 
uncertainty should be taken into account. For a 
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Salpcter (1955) slope, there is (roughly) a one-to-one 
correspondence between the halo normalization un- 
certainty and the corresponding predicted upper mass 
limit uncertainty. 

While we do agree with T97 that the lower limit 
to my is ~ 40 M© (or ~ 30 M Q , if we adopt the 
extreme A78 yields), our more realistic exploration 
of input physics "space" demonstrates that we sim- 
ply cannot constrain the upper limit to any useful 
accuracy. Taken together, we can only conclude that, 
by this technique alone, m\j > 40 M Q , for IMF 
slopes x = 1.3 — ► 1.6. Fixing the IMF slope to 
that of Salpcter (1955), we can only constrain m\j 
to lie somewhere between ~ 40 M Q and ~ 140 M Q . 
Finally, it would appear to be difficult to reconcile 
any my 100 M with either the A96zams or 
WW95;+0.60 halo normaliztion grids. 

While promising (provided existing discrepancies 
in Type II SNe yields arc eliminated), at the present 
time, unfortunately, this technique, by itself, does not 
substantially improve or constrain our understanding 
of the upper mass limit to the solar neighborhood 
IMF. 

We wish to thank Takuji Tsujimoto and Sylvia 
Becker for a number of helpful correspondences. BKG 
acknowledges the financial support of an NSERC 
Postdoctoral Fellowship. 
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